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ABSTRACT 

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the 
Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally 
developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find 
isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorpo- 
rates several advances, including a metric-based gridding of the search parameter space (frequency, frequency 
derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies 
between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs J 1803-2 149 
and J21 11-1-4606, are young and energetic Galactic -plane pulsars (spin-down power above 6x 10^^^ erg s"^ and 
ages below 100 kyr). The seven remaining pulsars, PSRs J0106-F4855, J0622-F3749, 11620^927, J1746-3239, 
J2028-F3332, J2030-F4415, J2139-F4716, are older and less energetic; two of them are located at higher Galac- 
tic latitudes {\b\ > 10°). PSR JO 106-1-4855 has the largest characteristic age (3 Myr) and the smallest surface 
magnetic field (2x lO^^G) of all LAT blind-search pulsars. PSR J2139-I-4716 has the lowest spin-down power 
(3x10^^ erg s~^) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency 
observations, only PSR JO 106-1-4855 has detectable pulsations in the radio band. The other eight pulsars belong 
to the increasing population of radio-quiet gamma-ray pulsars. 

Subject headings: gamma rays: stars, pulsars: general, pulsars: individual (PSR J0106-I-4855, 
PSR J0622-F3749, PSR J1620^927, PSR J1746-3239, PSR J1803-2149, 
PSR J2028-F3332, PSR J2030-F4415, PSR J21 11-^4606, PSR J2139-F4716) 
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1. INTRODUCTION 

The Fermi Gamma-ray Space Telescope has been operat- 
ing since launch in 2008 June. The Large Area Telescope 
(LAT) on board the Fermi satellite has an effective area of 
~ 0.8 (on- axis, above 1 GeV) and is sensitive to photons 
with e nergies from 20 MeV to more than 300 GeV fAtwoodl 
etal.|2009 ). Post-processing assigns an arrival time, energy E, 



we report on the discovery of nine pulsars from this ongoing 
effort in which a total of 109 sources selected from the Fermi- 



and direction to the photons; we call these "events". The ar- 
rival times are accurate to better than 1 /is, and the energies 
are accurate to within 15% between 0.1 and 10 GeV on-axis. 
The directional precision is energy-dependent: 68% of pho- 
tons have angular offset less than ~ 0.°8 x(£;/GeV)"*^-^ from 
the true direction ( Abdo et al. 2009c ). 

Gamma-ray pulsars are among the most interesting sources 
observed by the Fermi LAT. These are rapidly- spinning neu- 
tron stars whose regular beam of gamma-ray emissions passes 
by the detector with each rotation. The Fermi LAT has de- 
tected gamma-ray pulsations from more than 50 normal and 
millisecond pulsars (MSPs) previously discovered in radio- 
frequency searches (see e.g.; | Abdo et a l. 2010ci Ran som et al. 
2011} [Keith et aLpOlllfCognard et al.,2011.; Theureau et al. 



mi 



In contrast, so-called "blind" searches for gamma-ray pul- 
sars are not guided by any prior knowledge (e.g. from ra- 
dio or X-ray observations) of the pulsars' parameters. Pre- 
vious blind searches of the data recorded by the Fermi-LAI 
have been spectacul arly successful; they have discovered 26 
gamma-ray pulsars ( |Abdo et aL]|2009a[ |Saz Parkinson et aF 
'0101120111). This paper describes the blind" discovery of ^ 



2010[|2011| ). inis paper 
additional pulsars using a new method. 

The blind search problem is computationally demand- 
ing because the relevant pulsar parameters (typically sky- 
position, frequency /, and spin-down r ate /) are not known a 
priori and must be explicitly searched ( [Chandler et al.|200T[ ). 
So far, blind searches of LAT data hav e used a clever time - 
differencing t echnique" as described in|Atwood et al.[([2006|); 



Ziegler et"aLl ( [2008| . One powerful motivation for seeking 
even better methods is the application to MSPs, of which pre- 
vious blind searches have not found any. The blind search for 
MSPs is more computationally challenging due to the higher 
frequency range that must be covered. Moreover, most MSPs 
are in binary systems, where the orbital modulation parame- 
ters must also be searched, adding orders-of-magnitude to the 
computational complexity and challenge. 

This paper presents results from a new effort to find isolated 
gamma-ray pulsars (including MSPs), using a novel method 
inspired by computationally-efficient techniques recently de- 
veloped to search gravitational-wave detector data for weak 
continuous - wave signals from rap idly- spinning isolated neu- 
tron stars ( [Pletsch & Allen[[2009[ [PTetsch 2010^ ). In partic- 
ular, the search method uses a recently-aevelop ed optimal 
incoherent-combination method first described in ' Pletsch &I 
p.llen ( 200 9]) together with a "sliding coherence window tech- 
nique" (Pletsc h|20lT[ ). 

The method was originally intended to find isolated MSPs 
up to 1.4 kHz spin frequency with spin-down rates in the 
range -5x 10"^^ Hz s"^ < / < 0, with a corresponding range 
of characteristic ages r = -f/2f. However, the search tech- 
nique is also sensitive to normal (non-MSP) isolated pulsars, 
and several normal pulsars were discovered soon after we be- 
gan. This prompted us to extend the spin-down-rate search 
range, in order to include younger objects (having smaller r), 
down to ages (^ kyr) comparable to the Crab pulsar. Here 



LAT Second Source Catalog ( Abdo et aL][2011b| are being 
searched for new gamma-ray pulsars. 

The outline of this paper is as follows. Section |2| describes 
the LAT data preparation, and the selection of unidentified 
sources to search for previously unknown gamma-ray pulsars. 
The new hierarchical search method is explained in Section [3] 
and illustrated in Sectionpjwith a detailed example: our first 
discovery, PSR J 1620-4927. Section|5]presents the results for 
all the new gamma-ray pulsars. The search for counterparts in 
other regions of the electromagnetic spectrum is discussed in 
Section^ This is followed by a brief conclusion. 

2. SOURCE SELECTION AND DATA PREPARATION 



The Fermi-LAI Second Source Catalog (2FGL, Abdo et al. 



201 lb| ) lists 1873 sources, described by fits to elliptically 
shaped 95% -confidence sky regions. Among these sources 
576 are not associated with counterparts observed at other 
wavelengths and thus might contain unknown gamma-ray 
pulsars. 

In searching for new gamma-ray pulsars, it is important 
to identify and exclude sources that are blazars. Previously- 
obse rved Fermi-LXr gamma-ray pulsars (see e.g. [Abdo et al. 
[201 0c) have sharp cutoffs in their emission spectra at a few 
GeV, and stable gamma-ray fluxes. In contrast, blazars emit 
above 10 GeV, and their fluxes vary with time. An illustra- 
tion of the different spectral pr operties and variabil ity behav- 
ior can be seen in Figure 17 of [Abdo et al.[ ( [201 fb] ). Gamma- 
ray pulsars tend to have large curvature significances ("Sig- 
nif_Curve" in the 2FGL catalog, which gives the improve- 
ment in the quality of the spectral fit when changing from a 
power law to a curved spectral model) and small variability 
indices ("Variabihty_Index" in the 2FGL catalog, a measure 
of flux instability over time). In contrast, blazars tend to have 
small curvature significances and large variability indices. 

For the search for gamma-ray pulsations we select 2FGL 
sources with curvature- significance values greater than 4a 
(here and throughout the manuscript, a denotes the stan- 
dard deviation for a Gaussian distribution), and variabihty in- 
dices smaller than 41.6. We further select bright objects by 
choosing sources with detection significances ("Signif_Avg") 
greater than lOcr. Finally, we restrict the list to sources with 
no known associations ("unassociated sources") or associated 
with known supernova remnants (SNRs). By applying the 
above selection criteria to a preliminary version of the 2FGL 
catalog, we obtained a list of 109 2FGL sources to search 
for gamma-ray pulsations. (One source, 2FGL J062 1.9-1-3750 
was associated with an active galactic nucleus (AGN) in the 
final 2FGL catalog.) 

The Fermi Science Tools (STp] v9r23pl are employed to 
select the Fermi-LAI events for our search. Using gtselect, 
we take events from 2008 August 4 to 201 1 April 6, with re- 
constructed directions within 8° of the gamma-ray sources, 
energies above 100 MeV, and zenith angles < 100°. Only 
events belonging to the Pass 6 "Diffuse" class are retained, 
as those events have the highest probability of being photons 
( [Atwood et al.|2009) . For the event probability weighting de- 
scribed below we use the P6_V1 1 Instrument Response Func- 



See [http : / / f ermi . q sf c . nasa . qov/ssc/data/ 



analysis /scitools / overview . htm l| and |http : //fermiT 



gsf c . na sa . gov/ ssc/dat a/analysis /document at ion/ 



Cicerone / for details and tutorials. 
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tions (IRFs). Using the gtmktime tool, times when the rock- 
ing angle of the satellite exceeded 52° are excluded. We also 
require that DATA_QUAL and LAT_CONFIG are unity, and 
that the Earth's limb does not impinge upon the Region Of 
Interest (ROI). 

As d emonstrated in previous work ( |Bickel et al. [2008 [[Kerr 
201 la[ ), the sensitivity of gamma-ray pulsation searches can 
be improved by weighting the jih event with a probability 
G [0, 1] that it originated from a putative pulsar. The pho- 
ton weights Wj are calculated by using a full spectral model 
of the region around the gamma-ray source, and by exploiting 
the IRFs to provide background rejection which is superior to 
simple angular and energy data- selection cuts. For the first 
time in this paper such a photon-weighting scheme has been 
applied in a blind search; further details are given in Section[3] 

To calculate the weights Wj for each event from a given 
source, we perform likelihood spectral analyses using the 
pointlike tool (see Kerr 2011b , for a description). For each se- 
lected source, we construct a spectral model for the region by 
including all sources of the 2FGL catalog found within 8° of 
the selected source, using the spectral forms given in the cat- 
alog. The spectra of the selected sources are modeled as ex- 
ponentially cut-off power laws, typical of known gamma-ray 

pulsars, of the form Nq (E / GeV) exp (-E /Ec), where Nq is 
a normalization factor, F is the photon index and Ec is the cut- 
off energy. The Galactic and extragalactic diffuse emission 
and residual instrument background also enter the calcula- 
tion of the weights. The Galactic diffuse emission is modeled 
using the gll_iem_v02_P6_Vl 1 _DIFFUSE map cube, while 
the extragalactic diffuse and residual instrument backgrounds 
is modeled using the isotropic_iem_v02_P6_Vl 1 _DIFFUSE 
template (a detailed description of these backgroun d models 
can be found in Section 3 of Abdo et al.l pOlOa*)). These 
models are available for download at the Fermi Science Sup- 
port Centeip^ The tool gtsrcprob is then used to calculate 
the event weights Wj based on the best-fit spectral models ob- 
tained from the maximum likelihood analyses. 

3. THE NEW SEARCH METHOD 

In a year, the LAT detects of order 10^ photons from a typ- 
ical gamma-ray pulsar; in the same year, a typical pulsar ro- 
tates at least 10^ times around its axis. The blind-search prob- 
lem is to find a rotational-phase model $(0 = luift + ft^ /I) 
and a sky-position that match the Solar System Barycen- 
ter (SSB) arrival times t of the different photons, where $ 
denotes the rotation angle of the star about its axis, in ra- 
dians, measured from its starting position at ? = 0, and ob- 
served at the SSB. The signal hypothesis is that the photons' 
arrival times are "clustered" near specific "orientations" of the 
star (i.e., $(0 mod 27r deviates from uniformity on the inter- 
val [0,27r]). The null hypothesis is that the arrival times of 
the photons are a random Poisson process. In this paper we 
do not explicitly indicate the dependence of $ on /, / and 
sky position, but this dependence is important and implicit in 
many formulae below. 

To find a matching phase-model, a grid of "templates" in 
the four-dimensional parameter space of sky position and 
(/,/) is constructed. Note that the 2FGL catalog sky posi- 
tions of the targeted unassociated sources based on the spatial 
distribution of events are typically not precise enough for pul- 
sar searches. A search grid of sky points around this catalog 

http://fenni.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 



position is needed to reduce signal loss arising from imperfect 
correction of the Doppler shifts caused by the Earth's orbital 
motion around the SSB. The need for sky gridding is partic- 
ularly acute for MSP spin frequencies. Therefore, in contrast 
to previously published blind searche^ we grid a circular 
sky-region centered on the 2FGL catalog source location us- 
ing a radius which is 20% larger than the semi-major axis 
of the 95% confidence elliptical error region (given by the 
"Conf_95_SemiMajor" parameter). 

Unfortunately the number of templates (grid points) re- 
quired to discretely cover the entire four-dimensional search 
parameter space increases as a high power of the coherent 
integration-time. Hence a fully-coherent approach for several 
years of data is computationally impossible. Therefore, we 
employ a search strategy which is designed to achieve maxi- 
mum overall sensitivity at fixed computing cosj^ 

To efficiently scan through years of Fermi-LAI data for 
previously unknown gamma-ray pulsars, we use a so-called 
hierarchical search approach. This is analogous to hierarchi- 
cal methods used in searches for gravitational-wave pulsars 
([Schutz & Papa 1999; Papa et al. 2000; Brady & Creighton] 
[zOnPt lAbBotfet al.||2009a|bt [Cutler et al. 2005 ). In a first 
semi-coherent stage, we here adopt the optimal metric-based 
gridding methods described in Pletsch & Allen ( 2009 ) along 
with the sliding coherence window technique (Pletsch 2011 ). 
In a second stage, significant semi-coherent candidates are au- 
tomatically followed up in a fully-coherent analysis. Finally, a 
third stage further refines coherent pulsar candidates by using 
higher harmonics. Full details of the complete search scheme 
will b e presented in forthcoming work ( Pletsch & Guillemot[ 
20TT] ). 



Here we first describe the principle of the method, to firmly 
establish the analogy with the existing gravitational-wave lit- 
erature. Then we describe what is done in practice, which is 
mathematically equivalent (up to justifiable approximations) 
but computationally more efficient. 

In the first stage, a semi-coherent detection statistic S is 
computed for each template. We refer to S as "semi-coherent" 
because it is effectively the incoherent sum over several years, 
of terms which are coherent over several days. The coher- 
ent terms are the power in Fourier bins, calculated by treating 
each photon arrival as a delta function in time. 

Denoting the arrival time of the jih event (photon) at the 
SSB by tj, the coherent power in a (Gaussian) window 
centered at time r is defined by 



Pr 



E 

7=1 



(1) 



The sum is taken over all photons (here, N = 8000) in the data 

set; the effective window duration is / e~^^^ dr = T/ Vl. 
As described in Section[2] the weights Wj estimate the proba- 
bility that the photon comes from the selected source. 
To form the semi-coherent detection statistic S, the values 



'Zie gler et^ |2008\ argues correctly that for two-week data-stretches, 
sky gridding is not essential. However for data stretches of length compara- 
ble to a year or longer, sky gridding is necessary to avoid significant loss of 
signal-to-noise ratio. 

If this constraint is removed, then it is obvious that a more sensitive 
method exists. 
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of Pr are summed ("incoherently combined") 



2 r ^ 

^ 7=1 



W 



(2) 



Note that in this definition we have subtracted a constant 
(phase model indepdendent) term. Because it contains a 
Gaussian window, the integrand in Equation ^ falls off ex- 
ponentially at early and late time (large values of |t|). Thus 
the limits of integration can be taken as the entire real line 
T G (-00, oo); to good approximation this gives the same value 
as integrating only over the total observation interval (about 
975 days in this search). 

The semi-coherent detection statistic S is an incoherent sum 
of powers, which discards the phase information over time pe- 
riods longer than of order T. This uniform overlap maximizes 
the search sensitivity for fixed T and for fixed computational 
resources ( |Pletsch|2Qll| ). For computational efficiency, in this 
search we choose theN = 8000 photons with the highest prob- 
abilities (largest values of wj). 

To understand how to compute S efficiently, one can ex- 
plicitly evaluate Equation p]). Completing the squares in the 
product of the Gaussians ana carrying out the integration over 
T, one obtains 



(3) 



7=1 ^=1 



7=1 



Here the^ effective duration of the Gaussian window is 
J ^-7TT /T -j^ practice, to compute the semi-coherent 
power efficiently, we replace the Gaussian window in Equa- 
tion Q with a rectangular window of the same duration T, as 
givenbelow in Equation ([5]). In this search, the width of the 
rectangular window is T = 2^^ s (^ 6 days). 

The template grid in parameter space is the Cartesian prod- 
uct of a rectangular two-dimensional grid in (/, /) and a sky 
grid which has constant density when orthogonally projected 
onto the ecliptic plane. The problem of constructing efficient 
search grids has been intensively studied in the context of 
gravitational wave searches (see e.g. [B rady et al. 1998; Brady 
& Creigh ton]20QQl|Prix|2007[|Pletsch & Allen,2009, Pletsch 
^10 ) and we employ these concepts here. The values of fre- 
quency are equally spaced, separated by the FFT frequency- 
bin width A/ = 1 /T. In practice, to reduce the fractional loss 
in S for frequencies not coinciding with Fourier frequencies, 
we use a computati onally-efficient inte rpolation, referred to 
as "interbinning" ( Ransom et al.||2002] ). The spacing in the 
other three dimensions is determined by a metric which mea- 
sures the fractional loss in the expected value of S , that arises 
if the signal is not located exactly at a grid point (Balas ubra- 
[manian et al.^1996, Owen 1996; Prix 2007 ). I n spin-down /, 
we use a uniform grid spacing A/ = ^/llOm/n^T^. In this 
search m = 0.3 is the maximum tolerable fractional loss in 5", 



and from Pletsch & Allen (2009), 



7=1 



T2 



(4) 



where t = ^jtj/N is the mean photon arrival time. The de- 
scription of this grid can be found in [Pletsch & Allen| ( [2009] ) 
with a detailed derivation in |Pletsch (|2010). The grid in the 
sky is determined by the same metric, permitting a maxi- 



mum fractional loss m in the value of S. The spacing of the 
sky grid is determined by the Doppler- shift arising from the 
Earth's (more precisely, the Fermi satellite's) motion around 
the Sun. At the north Ecliptic pole the angular spacing is 
A^ = \/lmc/TifD, where c is the speed of light and £) is a 
baseline distance (defined below). When the entire sky grid 
is projected into the plane of the ecliptic, the grid-points are 
uniformly spaced on the plane ('Astone et aL [2002[ [Abbott 
[et al. 2009a b ). This angular spacing is similar to the an- 
gular spacing in the diffraction pattern of a two-slit system, 
where the wavelength is c/f and the separation of the two 
slits is the straight-line distance D between two points on the 
Earth's orbit about the Sun. If the coherent integration time T 
is less than half a year, then D = (998 s) c sin(7rr / 1 yr). If the 
coherent integration time T is greater than half a year, then 
D = (998 s) c is the diameter of the Earth's orbit about the Sun. 

To compute S efficiently, a time series is constructed and 
subsequently Fourier-transformed into the frequency domain. 
The time series contains T A /bw bins, where A/bw is the to- 
tal frequency bandwidth being searched at a time using com- 



plex heterodyning^jat the center of A/bw (see e.g. Patel et al 



(2010)). The time series is initialized to zero, then the values 
of WjWke~^^^ are added into the bins determined by the 



time differences Atjk = tj-tk, for all pairs of photons j,k for 
which < \Atjk\ < T; the bin index is obtained by round- 
ing the absolute value of the product A/bw ^tjk to its nearest 
integer value. 

Then the array is Fourier-transformed into the frequency 
domain (exploiting the FFT) to obtain S over the entire / grid. 
Up to an overall window-dependent normalization, one can 
write for given values of /, / and sky position, 

A^ 

S=Y, Qd^tjk I /T) Wj Wk (jH) , (5) 

7,^=1 

where the rectangular function Q(x) is unity if < x < 1 and 
vanishes otherwisj^ 

Although other aspects are different, the use of an FFT ap- 
plied to time differences is very similar to techniques pre- 
viously use d in blind searches of Fermi-LAI data ( Atwood] 
et al.|2006[ Equation 3 therein has a typo which is corrected 



m Equation 2 of [Ziegler et al.|| (2008 )). The method used in 
Atwood et al.l (I2006I) was the first application of this classic 



method (e.g. Bla ckman & Tukey|1958| ) to gamma-ray astron- 
omy (in estimating the power spectrum an approximate auto- 
covariance function is calculated using a maximum lag, and 
then Fourier- transformed). 

In contrast to previous searches, our method uses an optimal 
gridding of the parameter space for both the semi-coherent 
and coherent stages, as well as an automated follow-up, and 
incorporates the spin-down corrections in a way that permits 
heterodyning and highly-efficient code. 

The search was done on the 1680-node Atlas Computing 
Cluster ( jAulbert & Fehrm ann 2008) built around four-core 
processors with 8 GB of random-access memory; for these 
we used a heterodyning bandwidth A/bw = 256 Hz. Breaking 
the full frequency-range of the search into frequency-bands 
allows the computation to fit into memory, and also allows the 

Complex heterodyning is a procedure which shifts frequencies in time- 
series data by a fixed offset //, . This is accomplished by multiplying the time 
series by g'^^AO , shifting all frequencies by //,. 

Note that by symmetry S is real, because interchanging indices j and k 
is equivalent to complex conjugating the exponential factor. 
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Fig. 1. — Search results for the newly- discovered pulsar PSR 11620^927 (the large black dot in the lower- left of each panel). The left panel |(a)| shows the 
semi-coherent search results, representing about 2 CPU-years of computing on a single core. (Since the computing cost scales as the square of frequency, a search 
up to 64 Hz as in previous searches ( Abdo et al. 2009a Saz Parkinson et al. 2010 ) would have taken about 1.4 CPU-days on a single core.) The bottom-left panel 
shows the semi-coherent detection statistic 5 as a function of / and /, maximized over the sky-grid. The value of S is represented by the colorbar. A further 
maximization over / is shown to the right; a further maximization oyer / is shown above. PSR 11620^927, the darkest point near the bottom left of each figure, 
stands out clearly from the noise. In the same form, the right panel |(b) ( presents the fully-coherent follow-up search results of the previous candidate (and every 
other candidate "dot") shown in the left panel. The quantity plotted is now the fully-coherent detection statistic P over the entire data set. As explained in the 
text, for each candidate this covers a region of parameter space which is four steps of the semi-coherent grid in each dimension. 



use of different sky grids in each band. This further reduced 
the computational cost, since the number of required sky grid 
points increases with the square of frequency. 

After computing S on the four-dimensional grid in param- 
eter space, points with statistically-significant values of S are 
candidates for possible pulsar signals, and are followed up in 
a second stage. This is done by "refining the grid" and in- 
creasing the coherent integration time. This is an hierarchical 
scheme which is analogous to "zooming": successively swap- 
ping microscope objectives for ones of higher magnification, 
then re-ce ntering the interesting point on the slide (see e.g. 
Cutler et al. 2005 ; Krishnan et al. 2004). In our case, this is 



done by constructing the fully-coherent detection statistic P 
over the entire data set (or equivalently taking the Gaussian 
window-size T ^ oo in Equation ([T])) obtaining 



-mtj) 



(6) 



where for convenience we have normalized by the positive 
constant k given by 



1 ^ 



(7) 



7=1 



The computing cost to coherently follow up a single candidate 
is negligible in comparison to the cost of the previous semi- 
coherent search. 

In selecting statistically-significant semi-coherent candi- 
dates which are automatically followed up in the second stage 
using a fully-coherent analysis, we do not use a fixed thresh- 
old to define "statistical significance". In the semi-coherent 
stage, the search code keeps an internal list of the strongest 
signal candidates. Each member of this list is coherently fol- 
lowed up and corresponds to the largest value of S detected in 



eight adjacent spin-down values for the entire heterodyning 
frequency bandwidth and a single sky point. 

The refined grid of the fully-coherent follow-up covers a 
region of parameter space of size (4A/) x (4A/) x (4A6 x 
4A6) when projected into the ecliptic plane. In other words, 
it covers a region whose volume is 256 times larger than the 
volume of a fundamental cell in the original grid: its extent 
in each dimension of parameter space is four grid-intervals. 
The refined grid has a spacing given by the previous formulae 
for A/, A/, and A^, except that the coherence time T is set 
equal to the length of the entire data set, and 7=1. Since in 
this case only a small parameter- space region around the can- 
didate is explored, it is computationally-efficient to compute 
P directly in the time-domain (FFTs are not used), exploiting 
the sparsity of the photon data. 

If the value of P, which measures the fully-coherent power 
(in a single harmonic), is statistically significant, then in a 
third stage further refinement is carried out using higher har- 
monics (Fourier-components). We adopt the so-called //-test, 
which has been widely used in X-ray and gamma-ray pulsar 
detection (de Jager et al. 1989; de Jager & Biisching 2010]). 
This test measures the statistical significance of the energy m 
the first 20 (non-DC) Fourier-components of the pulse-profile 
as a function of phase. As in Equation ([T]), the H-tQst can 
also be m odified to include the photon probabilit y weig hts wj 
(see |Ken-|2011a| ). Note that Equation (5) in Ke^ pOll a) con- 
tains an error; corrected formulae used in this work are given 
below. 

The weighted //-test statistic is defined as follows. For each 
photon arrival time tj, the pulse phase xj (between zero and 
one) is calculated as xj = [^(tj) mod In] / In. The pulse pro- 
file (for <x < l)isa sum over the photons 



p(x) = ^Wj5(x-Xj), 



(8) 
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Fig. 2. — ■ The we ighte d //-test statistic for PSR J 1620-4927. Panels [(a)| and |(b)| show contour plots of the weighted //-test statistic as a function of sky 



position 



(a) and (/, /) |(b)| These are peaked at the parameter-s pace location indicated by the black cross and the axes show the offset from these values. Panel (c) 
shows how the maximum weighted //-test statistic [shown in |(a)| and |(b)) accumulates with time. The //-test increases (approximately) linearly with time, as 
expected for a pulsar that is emitting uniformly. 



where xj is the pulse phase of the jth photon and S(x) is a 
one-dimensional Dirac delta-function. It can be expressed as 
a Fourier series 



p(x) = i<i ^ aiQ 



(9) 



i=-oo 



which implies that the (complex) Fourier coefficients ai are 
given by 



1 



(10) 



7=1 



where the definition of k is identical to Equation ([Tf . The nor- 
malization of the Fourier coefficients has been chosen so that 
//"the photon arrival times are uniformly distributed, indepen- 
dent random variables, then in the limit of large numbers of 
photons, for £>0, ^(ai) and ^(a^) are independent Gaussian 
random variables with zero mean and unit variance. Finally, 
the weighted //-test statistic is defined as 



H = max 

1<L<20 



E 



4(L-1) 



(11) 



The quantity which is subt racted, 4(L- 1), is m otivated by an 
empirical numerical study ( de Jag er et al.|1989| ), providing the 
best omnibus test for unknown pulse profiles. 

Maximizing H over sky-position, frequency /, and spin- 
down rate / selects the "narrowest and sharpest" overall pulse 
profile. In contrast, maximizing the statistic P=\a\\^ of 
Equation ^ favors putting more power into the lowest har- 
monics. 

Using H as the test statistic, a further stage of parameter- 
space refinement is done in the same way as before: in the fre- 
quency, spin-down, and sky parameters, we cover regions that 
include four grid-steps (in each dimension) of the previous 
grid. The chosen refinement factor (ratio of the number of grid 
points after and before refinement, in each dimension) is about 
an order of magnitude. At each grid point, the weighted H- 
statistic of Equation ( fTT] ) is found. The parameter-space point 
with the largest statistic is selected as our best estimate of the 
pulsars' parameters which are then further refined through the 



timi ng-a nalysis procedure ( |Ray et al.|2011| ) described in Sec- 
tionl5?2l 

The hierarchical search pipeline has been validated by suc- 
cessfully recovering previously known gamma- ray pulsars, 
including some of the brightest gamma-ray MSPs ( |Abdo et 
[201 Oc; Pletsch & Guillemot 2011). In the next section, the 
complete search scheme is illustrated with a detailed exam- 
ple. 

4. EXAMPLE: RESULTS FOR PSR 11620^927 

Before giving the results for all of the pulsars that have 
been discovered with this new search method, we first go 
through a single example in detail. This illustration uses 
PSR 116 20^9 27, the first new pulsar found in this work. 

Figure [T(a)l shows the first stage of the analysis. For each 
/ and / value, the largest value of S found in the sky grid 
is displayed as an intensity. The point of highest intensity 
corresponds to PSR J 1 620-4927 's / and / parameters. 

The second analysis stage is th e automated follow-up of all 
candidates shown in Figure [T(a)l As explained in Section [3] 
this is accomplished by carrymg out a fully-coherent search 
over a small region of parameter space around each candi- 
date which has been identified as statistically significant in 
the previous semi-coherent stage. If the candidate found in 
the semi-coherent step was simply a statistical outlier, then 
the fullv-coherent statistic P will not be significant. Fig- 



ure | 1(b) [p resents the results of the fully-coherent statistic P for 
PSR 11620^927. Again the point of highest intensity com- 
pared to the background, indicating the presence of a coherent 
signal in the data set, is due to the new pulsar. 

The importance of the photon probability weights wj can 
be illustrated by repeating the analysis using the same 8000 
photons with the weights set to unity. (Note that the weights 
were used in selecting the 8000 photons, so this is n ot a 
complete comparison.) The result is that Figure [l(a)[ still 
shows a statistically-significant outlier in the semi-coherent 
search step, which is followed up automatically and gives a 
statisti cally- significant outlier in the fully-coherent output of 
Figure [l(b)[ In both steps the signal and its statistical sig- 
nificance are reduced, but the pulsar is still found. However 
for some of the other new pulsars, this is not the case: if the 
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weights are set to unity, then the pulsar is not detected. Over- 
all the weights play a larger role for sources in crowded re- 
gions of the sky, for example near the Galactic plane. 

At the third stage, further refinement is carried out by max- 
imizing the weighted //-test statistic over a small region of 
parameter space around the most significant candidate from 
the previous step. As described in Section [3] the parameter- 
space grid used at this stage is yet another oraer-of-magnitude 
finer than the one of the previous stage. Figure |2| shows the 
corresponding results for PSR 11620^927. 

The parameter- space point with the largest //-test statistic 
is selected as our best estimate of the pulsar parameters at this 
stage. The uncertainties in these estimated parameters can be 
obtained by using the fact that a Icr deviation for a Gaussian 
distribution has a value ^ 0.6 of the maximum. In the current 
example, as shown in Figure |2] the value of the maximum 
weighted //-te st statistic is 56a 

Figure |2(c)| illustrates the way that the weighted //-test 
statistic accumulates over the total observation time. The lin- 
earity of this plot indicates that the properties of the source, 
and of the data, were time- stationary; it shows that the pul- 
sar's emission is not changing with time. In this way it also 
gives additional confidence that the pulsar is real and not a sta- 
tistical noise outlier. (We have not provided analogous plots 
for the other eight pulsars reported in this paper, but they are 
similar.) 

As shown in de Jager & Buschingl ( [201Q| ) and |Kerr| ( [2Qlla| |, 
the weighted //-test has a false alarm probability approxi- 



mately described by PpA ~ ^ 



-0.4// 



For this example, the as- 



sociated false alarm probability is log^g/^FA ~ -98. We have 
not tried to rigorously estimate the trials factor, but since 
the total number of floating-point computations that could be 
performed by our computing systems over a period of some 
months is less than 10^\ we conclude that PSR 11620^927 
is a real gamma-ray pulsar and not a statistical outlier. 

5. RESULTS FOR ALL NEW PULSARS 

Table [T] shows the names and sky positions of the nine 
discovered pulsars, and also lists known source associations. 
The inferred rotational parameters of the new pulsars are pre- 
sented in Table |2| In addition, further derived parameters are 
given, including the spin-down power E = -An^Iff, where 
/ = 10"^^ g cm^ is the assumed neutron star moment of iner- 
tia, and the magnetic field strength at the neutron star sur- 



il9 



3 ,5. 



^)^/^ G and at the light cylinder 



facets ^3.2 X 10 

5lc ~ 2.94 X \Qi^{-ff s^)^/^ G, respectively 

Figure p] plots the newly discovered pulsars on a /-/ di- 
agram, wnere they can be compared with the known pulsar 
population. One can see that the new pulsars lie in a simi- 
lar region (and hence belong at large to the same population) 
as the previously-found blind-search Fermi-LAI gamma-ray 
pulsars. 

Two of the nine new systems, PSRs J 1803-2 149 
and J21 11-1-4606, are young energetic pulsars 
(^ > 6 X 10^^ erg s"^ and r < 100 kyr), located near the 
Galactic plane. Among the seven remaining less-energetic 
and older pulsars, five are also located near the Galactic plane 



< 5°) and two are found at higher Galactic latitudes 
> 10°). One of these, PSR 10106-^4855, has the largest 



characteristic age r = 3 Myr and the lowest surface magnetic 
field strength (5s ^ 2 x lO^^G) of afl LAT blind-search 
pulsars found to date. Also standing out, PSR J2139-I-4716 
has the smallest spin-down power (^ = 3 x 10^^ erg s"^) 
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Fig. 3. — The frequency / and spin-down rate / of the pulsar population. 
The black dot s show the approximately 1800 pulsars in the ATNF catalog 
(Manchester et al. 2005), excluding pulsars in globular clusters. The green 
dots show the 26 gamma-ray pulsars discovered in previous blind Ferm/-LAT 
searches (Abdo et al. 2009a Saz Parkinson et al.|2010t ,,2011). The nine red 
stars show the newly-discovered gamma-ray pulsars reported in this paper. 
The dotted lines indicate constant characteristic ages r, the dashed lines show 
contours of constant spin-down power E and the dashed-dotted lines signify 
contours of constant surface magnetic field strength B^. 



among all non-recyc led gamma-ray pulsars ever detected (cf . 
IDeLuca et al.pOllj ). 

5.1. Source associations 

As listed in Table [T] all but one of the new pulsars have 
gamma-ray coun terparts from the Fer mi-LAI First Source 
Catalog (IFGL, |Abdo et al.] |2010a| ). The exception is 
PSR J2028-F3332, which first appeared in the 2FGL catalog. 
In addition, these eight IFGL unassociated sources have been 
classified as likely pulsar candidates based strictly on a com- 
parison of their ga mma-ray propertie s with previously de- 
tected LAT pulsars ([Abdo et al.|2011a| ). 

PSR J2111-i-46(fe IS associated with the source 
OFGL J2110 8-F460 6 from the Fermi Bright Source List 
( [Abdo et d^ l2009b). Only ten objects from this catalog still 
remain to be identified. 

The new system PSR J 1620-4927 is associated w ith the 
unide ntified AGILE source lAGL 11624^946 ( Pittori et al.| 
20091. This AGILE source is however about ten times 
brighter, with a flux of (67 ± 13) x 10"^ ph cm"^ s"^ In ad- 
dition, it has an error radius of 0.°58 and encloses other 2FGL 
sources. The new gamma-ray pulsar could therefore only 
explain a fraction of the total flux detected by AGILE from 
lAGL 11624^946. 



While not formally associated in the 2FGL catalog (Abdo 
|eraLl|2QTTB , PSR 12028-^3332 lies within the 99% error 
contours for the EGRET source 3EG 12027-^3429. How- 
ever, the pulsar accounts for only about a fifth of the 
flux for the EGRET source. Another nearby LAT source, 
2FGLJ2025. 1-^3341, is associated with an AGN and has a 
peak LAT flux that can account for the remainder of the mea- 
sured EGRET flux. 

We have searched the TeVCat online catalog of TeV 
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source^ for very high-energy counterparts to the newly- 
discovered pulsars, but found no associations. This is not 
surprising, considering that most of the new gamma-ray pul- 
sars have moderate spin-down power (~ 10^"^ erg s"^) and rel- 
atively large characteristic ages (> 100 kyr), while pulsars 
associated with TeV Pulsar Wind Nebulae (PWNe) tend to 
be young and energetic. PSRs J1803-2149 and J21 11-^4606 
have properties that make associations with TeV PWNe plau- 
sible; however, no detections of TeV emission from these di- 
rections have been reported. 

5.2. Timing analysis 

For each new gamma-ray pulsar the definitive parameters as 
given in Tables H and are determined via precise measure- 
ments of pulse Times Of Arrival (TOAs) and by fitting the 
parameters of a timi ng model to thes e measurements using 
methods described in |Ray et aLlpOll >. 

To obtain a set of TOAs, photons are extracted for each 
source using a radius and minimum energy cut to optimize 
the signal-to-noise ratio for each pulsar. Then these data are 
subdivided into segments of about equal length. The set of 
best-guess pulsar parameters from the search (maximizing H) 
is used as a preliminary ephemeris to fold the photon ar- 
rival times, producing a set of pulse profiles. The TOAs are 
then measured by cross-correlating each pulse profile with a 
multi-Gaussian or kernel-density template derived from fit- 
ting the entire data set. This is done u sing the unbinne d max- 
imum likelifioodjnethoddes Then 
Tempo 2 i Hobbs et al.|2006| ) is used to fit the TOAs to a tim- 
ing model including sky position, frequency and frequency 
derivative. 

As often is the case, the youngest object found in this work 
(PSR J21 11-^4606, T = 17 kyr), exhibits particularly large 
timing noise. In turn, this requires including higher-order 
frequency-derivative terms to make a good timing model fit; 
in this case up to the fourth derivative of the frequency as 
shown in Table |2| For the same reason, the timing models in 
TableElfor PSRs J2030-F4415 and J1803-2149 include terms 
up to the third and second frequency derivative, respectively. 

Based on these timing solutions. Figures |4|-[l2] show the 
resulting phase-time diagrams and pulse promesior each of 
the newly-discovered pulsars. These plots are obtained from 
calculating the phase for each of the 8000 photons selected 
for the blind search and using the parameters listed in Ta- 
bles [T] and |2| and weighting each event by its probability of 
originating iTom the pulsar wj. The integrated pulse profiles 
(weighted pulse phase histograms) are constructed with a res- 
olution of 32 bins in phase per rotation. The la error bars in 
the integrated pulse profiles are statistical only and are given 
by where j runs over all events falling into the 

same phase bin. (Note that the formula for the fractional er- 
ror has the opposite sign of the exponential: -1/2.) 

As seen in Figure s[4|-[T2[ eight of the nine pulsars have pulse 
profiles with two peaks. For pulsars where the two peaks 
are separated by nearly one-half of a rotation it is possible 
to detect or discover the pulsar at the second harmonic (i.e. 
at twice the actual spin frequency of the pulsar). For pulsar 
candidates which were discovered with (apparently) single- 
peaked profiles, we tested for the true fundamental spin fre- 
quency by folding at the subharmonic of the discovery fre- 
quency. If the subharmonic is the correct frequency, the two 



resulting peaks may satisfy one or more of the following con- 
ditions: offsets that are measurably different than 0.5 in phase, 
significant differences in the integrated weighted counts un- 
der each peak, or significant differences in the shape of the 
peaks in different energy bands. With modest signal-to-noise 
ratios, such determinations are not always conclusive. For 
PSR JO 106-1-4855, our identification of the true period was 
subsequently confirmed by the detection of radio pulsations 
with a single peaked profile. For the one apparently single- 
peaked pulsar in our sample, PSR J2 139-1-47 16, none of the 
above tests yield strong evidence for the profile being double 
peaked at half the frequency. Additional data will be required 
to strengthen this conclusion. 

In order to further characterize each pulse profile, we fit the 
pulsars' weighted gamma-ray peaks to Lorentzian lines. The 
derived pulse shape parameters are listed in Table|3] Note that 
PSRs J0622-F3749 and J 1746-3239 show indications of sub- 
structures in their main gamma-ray peak. For these pulsars 
a single Lorentzian line function is used for fitting the main 
component. Apart from PSR JO 106-1-4855, which is detected 
at radio wavelengths (cf. Figure [T3f , there is no particular 
reference for absolute phase of these pulsars. For the other 
eight pulsars we have arbitrarily assigned the absolute phase 
reference such that the first peak occurs at a value of 0.1 in 
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phase. The gamma-ray pulse profiles shown in Figures 
along with the pulse-profile parameters listed in Table 
very similar to those of the previously-discovered gamma^ray 
pulsars in blind searches (Abdo et al. 2009a ; Saz Parkinson] 
[et al.|20 10; Abdo et al. 2010c), further supporting the theory 
that the gamma-ray emission consists of fan beams produced 
in the outer magnetosphere. 



5.3. Spectral parameters 

The spectral parameters for the new pulsars are obtained 
by fitting each phase-averaged spectrum with an exponen- 
tially cut-off power-law with a photon index F and a cutoff 
energy Ec. The results for each pulsar are listed in Table |4| In 
addition to F and which are explicit parameters of the fit. 
Table |4| also gives the important derived physical quantities of 
photon flux Fioo (in units of photons cm"^ s"^) and the energy 
flux Gioo (in units of erg cm"^ s"^) for events with energies 
between 100 MeV and 100 GeV 

Analogous to the pulse-profile properties, F and Ec mea- 
sured for the new pulsars are also similar to those observed for 
previously-detected gamma-ray pulsars (Abdo et al. 2010c|. 
This is not surprising, because (as described in Section|2|) tar- 
get sources for our search have been selected based on similar- 
ity of their spectral properties to known gamma-ray pulsars. 

The distance (3 kpc) for one of the new pulsars 
(PSR JO 106-1-4855) can be inferred based on the dispersion 
of the radio pulse measuring the free electron column density 
(see Section 6]T]for details). As the remaining eight pulsars 
are radio-quiet, this method cannot be used to estimate their 
distance. Furthermore, none of the pulsars is associated with 
a known supernova remnant, preventing us from deriving dis- 
tance estimates from such source associations. 

However, it is still possible to obtain a crude estimate of the 
distance to the new pulsars, by exploiting the observed corre- 
lation between the gamma-ray luminosity and the spin- 
down power E for other gamm a-ray pulsars with distance 
measures (cf. [Abdo et aLl|2010c| ). Based on this correlation 
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"pseudo gamma-ray luminosities" Lps are derived as 



Lps - 3.2 X 10^^(^/10'^ erg s-')'/^erg s 



34 



(12) 



where the E values are obtained from Table |2| Assuming a 
geometrical correct ion factor /r^ = 1 for the emission cone 
Walters et al.||200 9) for all gamma-ray pulsars, the rela- 
tion = Airfo^ GxQod^ is used to convert the energy flux and 
pseudo gamma-ray luminosity i nto a "pseudo distan ce" d^^. 



following Equation (2) of .Saz Parkinson et al.| ( |2Q10| ): 



1.6 



(^/ 10^4 erg s-i) 



1/4 



100 



/10"ii erg cm"2 s" 



1/2 



kpc. (13) 



For each of the new pulsars the resulting values for Lps and d^^ 
are shown in Table |4| Note that these estimated gamma-ray 
luminosities and distances are su bject to a number of caveats, 
detailed in Saz Parkin son et aL] ( [201 Op , and could differ sig- 
nificantly from the actual values. 

5.4. Why were the new pulsars not found in previous blind 
searches? 

To examine whether the nine pulsars found with this new 
method could be detected w ith previous methods, we apply 
the sa me search method (A twood et al.p006^ Ziegler et al.| 
2008| ) used to successfully discover t he 26 previously-found 
blind-search LA T gamma-ray pulsars ( [Abdo et al.|2009a[|Saz 



Parkinson et a l. 2010, 2011). We select input data as done m 
previous searches; events are selected based on a fixed ROI 



and m inimum-energy cut as described in Saz Parkinson et al. 



( |2010 >. We use the same coherence window size (r = 2^^s) as 



in the previous (and in the first stage of this paper's search ) 
search. No photon weights are computed or used. No sky- 
gridding is done in the first stage of the search: only the 
2FGL-catalog sky position is used. 

The previous blind search code recovers three of the nine 
pulsars (PSRs J1620-4927, J2028-F3332 and 12111-^4716) 
when the 2FGL sky locations are searched. Two of the re- 
maining six pulsars are detected only if the correct known sky 
position is used (as opposed to the 2FGL-catalog position). If 
in addition the ROI and energy cuts are optimized (scanning 
different values) then three of the remaining four pulsars are 
detected. Finally, if the coherence window size is doubled 
(which dramatically raises the computational burden in a full 
blind search) then the last pulsar is found. 

As compared to the previously-published blind-pulsar 
searches of LAT data, the methods used for this work in- 
corporate several significant improvements in sensitivity and 
computational efficiency, as well as a longer data set, that ex- 
plain why these new discoveries are made. First, the use of 
efficient parameter space gridding over both sky position and 
frequency derivative allows pulsars to be found that are much 
farther from the LAT catalog position than is possible with 
previous searches. In addition, using photon weights for both 
event selection and for the search computations ensures that 
the detection significance is near optimal with only a single 
trial over event selections. In contrast, methods that use a 
"cookie cutter" event selection must either search over two 
additional parameters (minimum energy and radius for the se- 
lection), or suffer a sensitivity penalty from potentially non- 
optimal event selections. A key factor is that while weighting 
the photons provides only a modest sensitivity improvement 
over optimal cookie cutter selections, the improvement can be 
large when compared to non-optimal cuts. 



One might think that the additional systems found in this 
paper come about purely from the significantly-increased 
computer power that was available. This is not the case: the 
improved methods deserve almost all of the credit. The im- 
proved methods are so computationally efficient that had we 
searched only up to 64 Hz as was done in previous searches, 
we could have searched all 109 selected 2FGL sources over 
the initial spin-down range (|/| < 5x10"^^ Hz s"^) and would 
have found seven of the nine new pulsars using less than 5000 
CPU hours. These are modest resources in comparison with 
those used in previous blind searcheq^ The remaining two 
pulsars would only have been found if the spin-down range 
were increased by a factor of twenty to |/| < 10"^^ Hz s"^ 
increasing the required CPU time to about 100000 hours. 

6. RADIO AND X-RAY COUNTERPART SEARCHES 

6.1. Pulsed radio emission 

These discoveries represent a substantial increase in the 
number of gamma-ray pulsars detected in blind searches of 
Fermi-LAI data. Of the 26 previously-discovered LAT blind- 
search puls^;s^_£nly3Jia^ 



band (ICamilo et al.||2009t |Abdo et al.1|2010bt [Saz Parkinson 



|et al.|2010|J201 lb and there are tigh t upper limits on the others 
dRay et al.l201l(lKeith et al.|2011|. In order to exploit these 
new discoveries in population studies of the relative beaming 
fraction and geometry of the radio and gamma-ray emission, 
it is essential to determine if they are also visible as radio pul- 
sars. 

Because the source list was chosen from Fermi-LAI pulsar- 
like unassociated sources, all of these sources have been 
previously searched for radio pulsations by t he Fermi Pul- 
sar S earch Consortium ( [Ransom et al.||20lT} [Hessels et al. 
201 1| ). We have re-analyzed all of these archival observa- 



tions, with increased sensitivity because we can now do a sin- 
gle frequency trial folding the radio data with the gamma-ray 
ephemeris, and search only in a single parameter, the disper- 
sion measure (DM). This greatly reduces the number of points 
searched in parameter space (the "trials factor") and implies 
that much smaller pulsed signal amplitudes are statistically 
significant. Where we saw an opportunity to go significantly 
deeper, we also made a number of new radio observations. 

The telescopes and observing configurations used are de- 
scribed in Table |5] and the characteristics of the individual 
radio observations are shown in Table |6l We compute the 
sensitivities using t he modified radi ometer equation given in 
Equation (A1.22) of [Lorimer & Kramer (2005): 



/3 



57; 



sys 



w 



Gy/n^tintAF V l/f- 



W 



(14) 



where (3 is the instrument-dependent factor due to digitiza- 
tion and other effects (when unknown, we assume /3 = 1.25), a 
value of 5 has been assumed for the threshold signal-to-noise 
ratio for a confident pulsar detection, Tsys = Tree + T^ky, G is 
the telescope gain, n^ is the number of polarizations used (2 
in all cases), tint is the integration time, AF is the observa- 
tion bandwidth, / is the pulsar spin frequency, and W is the 
pulse width (for uniformity, we assume W = 0.1/f). Note that 
conventionally the 3K background temperature is included in 
the receiver temperature Tree, which is measured on cold sky, 

32 The UCSC group has about 30000 CPU hours/month available; their 
previous search of 2 years of Fermi-LAI data used 65000 CPU hours. 
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and so Tsky represents the excess temperature from the Galac- 
tic synchrotron co mponent, which we estimate by scaling the 
408 MHz map of Haslam et al.| ( [1982 ) to the observing fre- 
quency V with a spectral index a = -2.6 (defined as *Si, oc z^"). 

We use a simple approximation of a telescope beam re- 
sponse to adjust the flux sensitivity in cases where the point- 
ing direction was offset from the true direction to the pulsar. 
This factor is given by ^ = ^-(^/hwhm) /i.5^ where 6 is the off- 
set from the beam center and HWHM is the beam half-width 
at half maximum. A computed flux density limit of *Smin at 
the beam center is thus corrected via division by q for targets 
offset from the pointing direction. 

For eight of the nine pulsars, we have established that they 
are indeed radio-quiet (or extremely radio-faint), as viewed 
from Earth. 

For PSR JO 106-1-4855 our observations and analysis have 
revealed very faint radio pulsations. In two of our archival 45- 
minute GBT observations at 820 MHz, we detect the pulsar at 
a DM of 70.87 ± 0.2 pc cm"^ with a single-trial significance 
of 4-5(7 in each observation. Although the detections are 
not individually very strong, we gain additional confidence in 
their veracity from the fact that the peak phases are consistent 
to 0.01 pulse periods when each observation is folded using 
the ephemeris determined from our gamma-ray timing. The 
summed radio profile shows a narrow peak with a duty cycle 
of ^ 2%, as shown in Figure 13 We estimate a flux density 
of ^ 20 /iJy in both observations, using the standard radiome- 
ter equation and a measurement of the off-pulse noise level. 
These radio pulsations are detected at a flux density below the 
nominal detection limit of 30 /iJy because the duty cycle is a 
factor of 5 smaller than the 10% duty cycle used in the sensi- 
tivity calculation for unknown pulse shapes. This corresponds 
to an equivalent flux density of 8 /iJy at 1400 MHz using a 
typical pulsar spectral index a = -1.7. As seen in Table |6l 
we have made two other observations of this source, one at 
350 MHz with the Green Bank Telescope and one at 1.4 GHz 
with Effelsberg. Neither of these observations detect the ra- 
dio pulsations. Accounting for this narrow duty cycle of the 
pulse, the minimum detectable flux density for the Effelsberg 
observation at 1.4 GHz was about 15 /iJy, so the non-detection 
in that observation is not surprising. This constrains the spec- 
tral index to be a < -0.5. On the other hand, a spectral index 
of -1.7 would imply a flux density at 350 MHz of 85 /iJy, 
which is above the nominal sensitivity of our 350 MHz ob- 
servation. This implies that either a > -1.3 or the sensitivity 
of that observation was affected by scatter or DM broadening 
or high er than expected sky background. Using the NE2001 
model ( [Cordes & Lazio| 20Q2 ), the measured DM corresponds 
to an estimated distance of 3.0 kpc, over a factor of 2 larger 
than the pseudo-distance. Given the radio detection, we can 
measure the so-called phase lag 5 = 0.073 ± 0.003 between 
the gamma-ray and radio emission. 

6.2. Continuum radio emission 

Pulsars and SNRs (Supernova Remnants) have the same 
origin, although the comparatively short lifetimes of SNRs 
means that the number of pulsar-SNR associations is quite 
small. However the rare associations are of high interest, as 
they constrain a number of pulsar and SNR parameters. A 
recent text-book example is the association of the g amma- 
ray pulsar J0007-F7303 w ith the shefl-type SNR CTA 1 ( [Abdo 
etal.|2008t|Sun etal.|20 11). 

For the nine new gamma-ray pulsars no association with 



a known SNR listed in the most recent SNR-catalog (|Green 



2009| ) is found. However, SNRs are difficult to identify m case 
they are faint, confused or distorted by interaction with dense 
clouds. We have used large-scale and Galactic -plane radio 
continuum surveys to search for structures in the vicinity of 
the new pulsars, which may indicate an association. 

For the two pulsars PSRs J0106-F4855 and J0622-F3749 lo- 
cated above 10° of Galactic latitude we used the MPIfR- 
survey sampleip^ t o extract maps from the 408 M Hz and 
1420 MHz surveys (Hasl am et al.|1982[[Reich|19 82), which 
show faint sources and extended smooth diffuse emission, but 
no discrete o bject. The high-resolution interferometric NVSS 
maps (Con don et al.|1998) show numerous compact sources, 
but no extended features. 

The remaining pulsars are located within 4° of Galactic lat- 
itude covered by Galactic-plane surveys. The southern- sky 
pulsar PSR 11620^927 is located towards the gradient of an 
extended em ission co mplex as seen i n southern-sky surveys 
dReich et al. 2001 ; Jo nas et al.||1998i). The South ern Galac- 
tic Plane Survey (SGPS, I Haverkorn et al.||2006] ) is insensi- 
tive to large scale emission, but shows no small scale struc- 
tures within 0.°5 of the pulsar. PSR J1746-3239 is located at 
b = -2?2 close to a emission ridge sticking out from the plane. 
The shell-type SNR G356.6-1.5 with a size of 20' x 15' ( |Gray| 
1 199 4) is part of this ridge and about 0.°6 away from the pul- 
sar. PSR J1803-2149 is located just 2' away from the peak 
of a flat spectr um 6 Jy therma l source visible in all Galactic- 
plane surveys. [ Quireza et al.| #2006 ) list details for the HII- 
region with a deconvolved Gaussian size of 2^7 and a dis- 
tance ofl4 kpc. Also CO-emission was observed (S coville] 
et al. II 1987]). The pulsar is located at the periphery of this 



emission complex, but further studies are needed to investi- 
gate this region. PSR J2028-I-3332 is located at the southern 
boundary of the thermal Cygnus-X complex. The 1 1 cm sur- 
vey (Reich et al. 1984) shows patchy structures in the field 
and an extragalactic 1.8 Jy source 40' distant from the pul- 
sar. PSR J2030-I-4415 is seen towards the northern periphery 
of strong complex emission from the Cygnus-X region. Ded- 
icated studies are needed to find any emission associated with 
the pulsar. PSRs J21 11-^4606 and J2139-F4716 are both lo- 
cated in low emission areas, where the 21 cm and 11 cm Ef- 
felsberg Galactic-plane surveys show faint structures close to 
the noise level. 

6.3. X-ray 

Subsequent to the pulsar discoveries, we searched for 
archival X-ray observations covering the new pulsars' sky 
locations. As listed in Table [7] we have found short 
Swift observations (3-10 ks exposure) for five of the pul- 
sars. In addition there is a 6 ks-long XMM -Newton obser- 
vation for PSR J1620-4927 and a 23 ks-long Suzaku obser- 
vation for PSR JO 106-1-4855. For the two remaining pulsars, 
PSRs J1746-3239 and 12028-^3332, 10 ks-long Swift obser- 
vations were carried out following the discoveries. 

These X-ray data were analyzed, and no X-ray counterparts 
were found for any of the new gamma-ray pulsars. 

In order to estimate upper limits on the X-ray flux 
for each new pulsar, a power-law spectrum with a pho- 
ton index of 2 and a signal-to-noise of 3 is used. The 
values for the absorbing columns are estimated analo- 



gously to IMarelli et aT] ( |20T11 . In the 0.3-10 keV en- 
ergy range, the derived upper limits on the X-ray flux of 



^^ jhttp : //www.mpif r . de/ survey . html 
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the pulsars are between 1 and 3x10"^^ erg cm"^s"^ except 
for PSR 11620^927, where an upper-hmit X-ray flux of 
6x 10-14 gj-g cm-^s-i is obtained, and PSR J0106+4855, with 
an upper-limit X-ray flux of 2.3x10-1"^ erg cm-^s-^ These 
results, listed in Table [7] are consistent with the gamma- 
to-X-ray flux-ratios for previously-found Fermi-LAI pulsars 
( Marefli etall20TT] ). 

7. CONCLUSION 

This work reports on the discovery of nine gamma-ray pul- 
sars through the application of a new blind-search method to 
about 975 days of Fermi-LAI data. The new pulsars were 
found by searching unassociated sources with typical pulsar- 
like properties selected from the Fermi-LAI Second Source 
Catalog ( Abdo et al |2011bf . 

The sensitivity of blind searches for previously unknown 
gamma-ray pulsars is limited by finite computational re- 
sources. Thus efficient search strategies are necessary max- 
imizing the overall search sensitivity at fixed computing 
cost. We have developed a novel hierarchical sear ch method, 
adapted f rom an optimal semi-coherent method ( |Pletsch & 
[Alien 2009 ) together with a sliding coherence window tech- 
nique (Pletsch 2011) originally developed for detection of 
continuous gravitational-wave signals from rapidly- spinning 
isolated neutron stars. The first stage of the method is semi- 
coherent, because coherent power computed using a window 
of 6 days is incoherently combined by sliding the window 
over the entire data set. In a subsequent stage, significant 
semi-coherent candidates are automatically followed up via 
a more sensitive fully-coherent analysis. The method ex- 



tends t he pioneering methods first described in Atwood et al. 
( [20061 ). 

The new method is designed to find isolated pulsars up 
to kHz spin frequencies, by scanning a template grid in the 
four-dimensional parameter space of frequency, frequency- 
derivative and sky location. The optimal and adaptive sky 
gridding (not done in previously-published blind searches) is 
necessary, particularly at the higher search frequencies, be- 
cause the source-catalog sky -positions are not precise enough 
to retain most of the signal-to-noise ratio in year-long data 
sets. A fundamental new element of the method is the 
exploitation of a parameter-space metric (well- stu died in 
the continuous gravitational-wa ve context (see e.g. [Brady 
& Creighton][2^ [Pnxl 12007) )) to build an efficient tem 



plate grid, as well as a metric approa ch to construct an op- 
timal semi-coherent combination step ( [Pletsch & All en 20W ; 
[Pletsch 20 10 ). A further enhancement over previous searches 
is the sub-division of the total search frequency range into 
bands via complex heterodyning. This accommodates mem- 
ory limitations, parallelizes the computational work, and per- 
mits the use of efficient sky grids adapted to the highest fre- 
quency searched in each band. A photon probability weight- 
ing scheme (Kerr 2011a) is also used for the first time in a 
published blind search, further improving the search sensitiv- 
ity. 

The nine newly-discovered pulsars increase the previously- 
known Fermi-LAI blind-search pulsar population ( [Abdo et aT 
2009a[ [Saz Parkinson e t al.||2010| |2011| ) by more than one- 
third, aM brings the total number to 35. The inferred pa- 
rameters of the new pulsars suggest that they belong to the 
same general population as the previously-found blind-search 
gamma-ray pulsars. Deep follow-up observations with ra- 
dio telescopes have been conducted for all of the new pul- 
sars, but significant radio pulsations have only been found for 



PSR JO 106-^4855. The null results for the other eight pulsars 
indicate that they belong to the growing population of radio- 
quiet gamma-ray pulsars, which can only be detected via their 
gamma-ray pulsations. 

The computational work of the search has been done on 
the 6720-CPU-c ore Atlas Computing Cluster (Aul bert &^ 
Fehrmann|2008[ ) at the Albert Einstein Institute in Hannover. 
Recently, in August 2011, we moved the computational bur- 
den of the search onto the volunteer distributed computing 
system Ein stein® Hom^ ("Abbott et al.[ [2009a|b[ [Knispel' 
et al. 12010] ). This will provide significantly more computing 
power, and will allow a complete search of the parameter- 
space up to kHz pulsar spin frequencies. We also hope that in 
the future an improved version of these methods can be used 
to carry out blind searches for gamma-ray pulsars in binary 
systems. 

The combination of improved search techniques and much 
more powerful computational resources leave us optimistic 
that we can find still more gamma-ray pulsars in the Fermi- 
LAI data. These advances should also greatly increase the 
chance of finding the first radio-quiet gamma-ray MSP with 
the Fermi-LAI. We hope that further discoveries, and further 
study of these systems, will eventually provide important ad- 
vances in our understanding of pulsars, and of their emission 
mechanisms and geometry. 
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TABLE 1 

Names and Sky Locations of the Discovered Gamma-Ray Pulsars 



Pulsar Name 


Source Association 


R.A.^ 


Decl.^ 










(hh:mm:ss.s) 


(dd:mm:ss.s) 


(deg) 


(deg) 


T0106+485S 

jyj ± yjyj i i o ^ / ^ / 


2FGL J0106. 5+4854 
1FGLJ0106.7+4853 


01 ■06-25 06ri^ 


+48-55-51 8r2^ 


125.5 


-13.9 


J0622+3749 


2FGLJ062 1.9+3750 
IFGT T0622 2+3751 


06:22:10.51(2) 


+37:49:13.6(9) 


175.9 


11.0 


11620^927 


2FGLJ1620.8-4928 
1FGLJ1620.8-4928C 
lAGL 11624^946 


16:20:41.52(1) 


^9:27:37.1(3) 


333.9 


0.4 


J1746-3239 


2FGLJ1746.5-3238 
1FGLJ1746.7-3233 


17:46:54.947(8) 


-32:39:55.8(7) 


357.0 


-2.2 


J1803-2149 


2FGLJ1803.3-2148 

1FGLJ1803.1-2147C 

1AGLJ1805-2143 


18:03:09.632(9) 


-21:49:13(4) 


8.1 


0.2 


J2028+3332 


2FGLJ2028.3+3332 
3EGJ2027+3429 


20:28:19.860(5) 


+33:32:04.36(7) 


73.4 


-3.0 


J2030+4415 


2FGLJ2030.7+4417 
1FGLJ2030.9+4411 


20:30:51.35(4) 


+44:15:38.1(4) 


82.3 


2.9 


J21 11+4606 


2FGLJ21 11.3+4605 
1FGLJ21 11.3+4607 
0FGLJ21 10.8+4608 


21:11:24.13(3) 


+46:06:31.3(3) 


88.3 


-1.4 


J2139+4716 


2FGLJ2139.8+4714 
1FGLJ2139.9+4715 


21:39:55.95(9) 


+47:16:13(1) 


92.6 


-4.0 



Note. — A list of the nine new pulsars reported in this work, showing their sky locations 
and associations with cataloged gamma-ray source s. The associatio ns listed include sources from 
the Fermi-LKT Second Source Catalog (2FGL, Abdo et al |2011b|, the Fermi-LAI First S ource 
Catalog (IFGL, [Abd o et al. 20 10a), the FermPL AT Bright Source List (OFGL, Abdo eTaL] 
2009bt, the AG/LFT atalog (lAGEr |Pittori et al.t2 009 i, and the Third EGRET Catalog (3EG7 
kartman et al.TTW ). 

^ Right ascension (J2000.0) and declination (J2000.0) obtained from the timing model, where the 

numbers in parentheses are the statistical Icr errors in the last digits. 

^ Galactic longitude (/) and latitude (b), rounded to the nearest tenth of a degree. 



TABLE 2 

Measured and Derived Parameters of the Discovered Gamma-Ray Pulsars 



Pulsar Name 


/ 


/ 


Weighted 


r 


E 


Bs 


Blc 




(Hz) 


(-10-13 Hz s-i) 


//-test 


(kyr) 


(10^4 erg s-i) 


(10^2 G) 


(kG) 


J0106+4855 


12.02540173638(8) 


0.61881(7) 


843.1 


3081.1 


2.9 


0.2 


3.0 


J0622+3749 


3.00112633651(5) 


2.28985(4) 


288.8 


207.8 


2.7 


2.9 


0.7 


11620^927 


5.81616320951(5) 


3.54782(4) 


566.4 


259.9 


8.1 


1.4 


2.4 


J1746-3239 


5.01149235750(3) 


1.64778(3) 


249.8 


482.2 


3.3 


1.2 


1.3 


J 1803-2 149 


9.4044983174(2) 


17.25894(6) 


451.9 


86.4 


64.1 


1.5 


11.0 


J2028+3332 


5.65907208453(2) 


1.55563(2) 


1108.3 


576.8 


3.5 


0.9 


1.5 


J2030+4415 


4.4039248637(5) 


1.2576(2) 


584.8 


555.2 


2.2 


1.2 


1.0 


J21 11+4606 


6.3359340865(4) 


57.4218(3) 


554.3 


17.5 


143.6 


4.8 


11.1 


J2 139+47 16 


3.5354509962(2) 


0.2232(2) 


351.1 


2511.5 


0.3 


0.7 


0.3 



Note. — The reference epoch for all measured rotational parameters is MJD 55225 and the time range for all 
timing models is MJD 54682 - 55719. The derived quantities in columns 5-8 are based on the / and / values 
obtained from the timing model and are rounded to the nearest significant digit. To model the timing noise present 
in PSR Jl 803-2 149, a second frequency derivative is necessary: / =7.3(8) x 10-^^ Hz s-^. To model the timing noise 
present in PSR J2030+4415, higher frequency derivatives up to third order are necessary: / = -1.5(3) x 10-^^ Hz s-^ 
and / = -6(2) X 10-31 Hz s~^. To model the timing noise present in PSR J21 11+4606, higher frequency derivatives 
up to fourth order are necessary: /=2.30(5) x 10-^2 Hz s-^, / = -7.9(2) x 10-^0 Hz s-^ and '/ =3.2(4) x 10-^^ Hz s-^. 
The numbers in parentheses are the statistical la errors in the last digits. 
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JO 106+4855 




Pulse phase Pulse phase 



Fig. 4. — Phase-time diagram and pulse profile for PSR J0106+4855. The 
left panel shows the pulse phase at the arrival time of each photon, where 
the gray-scale intensity represents the photon probability weight. The upper 
right plot shows the summed probability weights: the integrated pulse profile 
using a resolution of 32 bins per rotation. The error bars represent the la 
statistical uncertainties. The four plots below resolve the integrated pulse 
profile according to separate energy ranges. For clarity, the horizontal axis 
shows two pulsar rotations in each diagram. In obtaining these plots the 8000 
events with the highest probability weights have been used (as in the blind 
search). 



J0622+3749 




Fig. 5.— Phase-time diagram and pulse profile for PSR J0622+3749. The 
plots have identical form as those shown in Figure |4] 



J1620-4927 




Fig. 6.— Phase-time diagram and pulse profile for PSR J 1620-4927. The 
plots have identical form as those shown in Figure l4l Note that for this par- 
ticular pulsar, based on the photon probability weights, no selected events 
(among the 8000) have energies in the 0.1-0.3 GeV range. 



J1746-3239 




Fig. 7.— Phase-time diagram and pulse profile for PSR J1746-3239. The 
plots have identical form as those shown in Figure[4] 
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Fig. 8.— Phase-time diagram and pulse profile for PSR J1803-2149. The 
plots have identical form as those shown in Figure |4] 



Fig. 1 0.— Phase-time diagram and pulse profile for PSR J2030+4415. The 
plots have identical form as those shown in Figureji] 



J2028+3332 




Pulse phase Pulse phase 



J21 11+4606 




1 L.J Z 

Pulse phase Pulse phase 



Fig. 9.— Phase-time diagram and pulse profile for PSR J2028+3332. The 
plots have identical form as those shown in Figure |4l 



Fig. 11.— Phase-time diagram and pulse profile for PSR J21 11+4606. The 
plots have identical form as those shown in Figurel4] 
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Fig. 12.— Phase-time diagram and pulse profile for PSRJ2139+4716. The 
plots have identical form as those shown in Figure |4] 
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TABLE 3 

Pulse-Profile Parameters of the Discovered Gamma-Ray Pulsars 



Pulsar Name 


Peak Multiplicity 


FWHMi 


FWHM2 




A 




J0106+4855 


2 


0.02 


± 


0.01 


0.05 ± 


0.01 


0.50 


± 


0.01 


J0622+3749 


2 


0.21 


± 


0.04 


0.03 ± 


0.03 


0.47 


± 


0.03 


J 1620-4927 


2 


0.13 


± 


0.04 


0.15 ± 


0.05 


0.21 


± 


0.03 


J1746-3239 


2 


0.32 


± 


0.07 


0.11 ± 


0.11 


0.63 


± 


0.06 


J1803-2149 


2 


0.04 


± 


0.02 


0.05 ± 


0.03 


0.40 


± 


0.02 


J2028+3332 


2 


0.10 


± 


0.02 


0.04 ± 


0.01 


0.38 


± 


0.02 


J2030+4415 


2 


0.09 


± 


0.03 


0.04 ± 


0.01 


0.49 


± 


0.02 


J21 11+4606 


2 


0.08 


± 


0.02 


0.07 ± 


0.02 


0.31 


± 


0.02 


J2139+4716 


1 


0.13 


± 


0.03 













Note. — For each of the nine pulsars, we give the parameters describing 
the shape of the pulse profile, including the peak multiplicity, the Full- Widths at 
Half Maxima (FWHM) of the peaks, and the separation A between the gamma-ray 
peaks for pulsars with more than one peak. 



TABLE 4 

Spectral Parameters of the Discovered Gamma-Ray Pulsars 



Pulsar Name 


r 


Ec 








dps 






(GeV) 


(10 ^ photons cm ^ s 


^) (10-11 erg cm-2 s-i) 


(10^3 ergs-i) 


(kpc) 


J0106+4855 


1.47 ±0.23 ±0.12 


3.31 ±0.92 ±0.08 


2.56 ± 0.77 ± 0.32 


2.40 ±0.31 ±0.04 


5.5 


1.4^^ 


J0622+3749 


0.59 ± 0.34 ± 0.09 


0.60 ±0.13 ±0.04 


2.21 ± 0.35 ± 0.08 


1.69 ±0.15 ±0.05 


5.3 


1.6 


J 1620-4927 


1.01 ±0.18 ±0.05 


2.44 ± 0.42 ± 0.28 


9.61 ± 1.68 ±0.94 


13.5 ± 1.0 ± 1.7 


9.1 


0.7 


J1746-3239 


1.33 ±0.08 ±0.35 


1.65 ±0.12 ±0.52 


9.97 ± 0.94 ± 1.85 


7.86 ± 0.41 ± 0.77 


5.8 


0.8 


J1803-2149 


1.96 ±0.11 ±0.20 


5.73 ± 1.72 ±2.07 


20.7 ±3.1 ±0.5 


13.1 ± 1.1 ± 2.1 


25.6 


1.3 


J2028+3332 


0.86 ± 0.21 ± 0.07 


1.53 ±0.24 ±0.08 


5.12 ±0.87 ±0.44 


6.09 ±0.41 ±0.13 


6.0 


0.9 


J2030+4415 


1.89 ±0.14 ±0.22 


2.16 ± 0.65 ± 0.67 


13.3 ± 1.4 ±0.2 


7.06 ± 0.48 ± 0.66 


4.8 


0.7 


J21 11+4606 


1.63 ±0.14 ±0.05 


5.43 ± 1.80 ± 1.56 


4.39 ± 0.69 ± 0.02 


4.13 ±0.34 ±0.30 


38.4 


2.7 


J2139+4716 


0.80 ± 0.27 ± 0.02 


1.02 ±0.21 ±0.07 


2.65 ±0.44 ±0.19 


2.51 ±0.21 ±0.01 


1.8 


0.8 


Note. — This table describes the spectral properties of each of the nine pulsars. 


modeling each spectrum as 


an exponentially cut-off 



power-law with photon indices F and cutoff energies Ec. The spectral parameters listed here for each pulsar are obtained from maxi- 
mum likelihood fits. The first quoted uncertainties are statistical, while the second are systematic and correspond to the differences in the 
best-fit parameters observed when doing the spectral analyses with the P6_V3 IRFs and associated diffuse emission models (namely, the 
gll_iem_v02 map cube and isotropic _iem_v02 template). For each object, the pseudo gamma-ray luminosity Lps and the pseudo distance dps 
are inferred from the apparent spin-down power E and the energy flu x Giqq above 100 MeV. Not e that these estimated gamma-ray lumi- 
nosities and distances are subject to a number of caveats, detailed in \Saz Parkinson et al.\\20ld\ , and could differ significantly from the 
actual values. 

^ Photon flux measured above 100 MeV. 
^ Energy flux measured above 100 MeV. 

^ The actual distance is 3.0 kpc, as inferred from the dispersion of the radio pulse measuring the free electron column density; see Section[6!T] 



TABLE 5 

Definition of Radio Observing Codes 



Obs Code 


Telescope 


Gain 


Frequency 


Bandwidth AF 


13- 


np 


HWHM 


Tree 






(K/Jy) 


(MHz) 


(MHz) 






(arcmin) 


(K) 


GBT-350 


GBT 


2.0 


350 


100 


1.05 


2 


18.5 


46 


GBT-820 


GBT 


2.0 


820 


200 


1.05 


2 


7.9 


29 


GBT-S 


GBT 


1.9 


2000 


7001^ 


1.05 


2 


3.1 


22 


Eff-Ll 


Effelsberg 


1.5 


1400 


250 


1.05 


2 


9.1 


22 


Eff-L2 


Effelsberg 


1.5 


1400 


140 


1.05 


2 


9.1 


22 


Jodrell 


Lovell 


0.9 


1520 


200 


1.05 


2 


6.0 


24 


AO-327 


Arecibo 


11 


327 


25 


1.12 


2 


6.3 


116 


AO-Lwide 


Arecibo 


10 


1510 


300 


1.12 


2 


1.5 


27 


Parkes-BPSR 


Parke s 


0.735 


1352 


340 


1.05 


2 


7.0 


25 



Note. — The sky locations of all nine pulsars have been searched for pulsating radio 
emissions. This Table gives the radio telescope and back-end parameters used in those 
observations, which are described in Table|6| 
- Instrument-dependent sensitivity degradation factor. 

^ The instrument records 800 MHz of bandwidth, but to account for a notch filter for RFI 
and the lower sensitivity near the band edges, we use an effective bandwidth of 700 MHz 
for the sensitivity calculations. 
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TABLE 6 

Radio Search Observations of the New Gamma-Ray Pulsars 



Target 


Obs Code 


Date 


'int 


R.A.^ 


Decl.^ 


Offset 


T 

^ sky 


c . 

'--^min 








(min) 


(J2000) 


(J2000) 


(arcmin) 


(K) 




J0106+4855 


GBT-350 


2009-10-25 


32 


01:06:37.7 


48:54:11 


2.7 


49.3 


136 




GBT-820 


2010-11-17 


45 


01:06:35.5 


48:55:30 


1.8 


5.4 


30^ 




GBT-820 


2010-12-17 


45 


01:06:35.5 


48:55:30 


1.8 


5.4 


30^ 




Eff-L2 


2011-06-01 


45 


01:06:25.1 


48:55:52 


0.0 


1.3 


31 


J0622+3749 


GBT-350 


2009-10-27 


32 


06:22:05.5 


37:51:07 


2.2 


46.2 


131 




Eff-Ll 


2010-05-14 


32 


06:22:14.7 


37:51:49 


2.8 


1.3 


30 




Eff-Ll 


2010-02-06 


10 


06:22:15.0 


37:51:48 


2.8 


1.3 


53 




Eff-Ll 


2010-07-10 


52 


06:22:13.0 


37:50:36 


1.5 


1.3 


22 




Eff-Ll 


2010-07-10 


55 


06:22:13.0 


37:50:36 


1.5 


1.3 


22 




GBT-820 


2010-12-12 


45 


06:21:59.0 


37:51:36 


3.3 


5.0 


32 




GBT-820 


2010-12-17 


45 


06:21:59.0 


37:51:36 


3.3 


5.0 


32 


J 1620-4927 


Parkes-BPSR 


2009-08-03 


270 


16:21:05.5 


-49:30:32 


4.9 


16.9 


42 




Parkes-BPSR 


2010-11-18 


144 


16:20:43.5 


-49:28:24 


0.9 


16.9 


42 




Parkes-BPSR 


2011-05-10 


72 


16:20:41.3 


-49:27:36 


0.0 


16.9 


58 


J1746-3239 


GBT-S 


2009-12-23 


60 


17:46:47.9 


-32:36:22 


3.8 


5.1 


30 




GBT-820 


2010-11-14 


45 


17:46:41.0 


-32:36:18 


4.6 


51.6 


85 




Parkes-BPSR 


2011-05-10 


72 


17:46:54.9 


-32:39:55 


0.0 


14.1 


54 


J1803-2149 


Eff-Ll 


2010-02-13 


25 


18:03:12.0 


-21:47:27 


1.6 


17.8 


55 




Eff-Ll 


2010-05-22 


32 


18:03:11.7 


-21:47:28 


1.5 


17.8 


48 




GBT-S 


2010-09-04 


65 


18:03:11.7 


-21:47:28 


1.5 


7.1 


13 


J2028+3332 


GBT-820 


2009-08-13 


60 


20:27:48.0 


33:32:24 


6.6 


10.3 


46 




GBT-S 


2010-09-20 


30 


20:28:18.0 


33:33:23 


1.4 


1.0 


15 




GBT-820 


2010-11-22 


45 


20:28:19.0 


33:32:53 


0.8 


10.3 


33 




GBT-820 


2010-12-17 


45 


20:28:19.0 


33:32:53 


0.8 


10.3 


33 




AO-Lwide 


2011-05-21 


45 


20:28:19.9 


33:32:06 


0.0 


2.1 


4 




AO-327 


2011-05-30 


25 


20:28:19.9 


33:32:06 


0.0 


112.9 


142 


J2030+4415 


Eff-Ll 


2010-02-07 


10 


20:30:55.0 


44:11:52 


3.8 


4.0 


62 




Eff-Ll 


2010-05-15 


32 


20:30:55.3 


44:11:53 


3.8 


4.0 


35 




Eff-Ll 


2010-07-10 


60 


20:30:59.2 


44:15:33 


1.4 


4.0 


23 




Eff-Ll 


2010-07-30 


60 


20:30:59.2 


44:15:33 


1.4 


4.0 


23 




GBT-820 


2010-11-22 


45 


20:30:54.7 


44:16:08 


0.8 


16.0 


38 




GBT-820 


2011-05-28 


183 


20:30:51.3 


44:15:38 


0.0 


16.0 


19 




Eff-L2 


2011-06-01 


45 


20:30:51.5 


44:15:37 


0.0 


4.0 


35 


J21 11+4606 


GBT-820 


2009-09-19 


60 


21:11:22.8 


46:05:53 


0.3 


16.0 


33 




Jodrell 


2011-06-22^^ 


14x60 


21:11:24.0 


46:06:29 


0.0 


9.0 


14 


J2139+4716 


Eff-Ll 


2010-07-10 


60 


21:39:52.3 


47:13:43 


2.6 


2.0 


22 




GBT-350 


2009-10-25 


32 


21:39:53.2 


47:15:22 


1.0 


74.9 


171 




GBT-820 


2010-12-11 


45 


21:39:53.5 


47:13:30 


2.7 


8.2 


34 




GBT-820 


2010-12-18 


45 


21:39:53.5 


47:13:30 


2.7 


8.2 


34 




Eff-Ll 


2010-05-15 


32 


21:39:56.9 


47:15:28 


0.8 


2.0 


29 



Note. — The sky locations of all nine pulsars have been searched for radio pulsations. Only for 
PSR J0106+4855 radio pulsations are detected. The minimum detectable flux density <Smin for each 
observation is computed at the observing frequency using Equation ([14} and the parameters in Table 
|5] as described in the text. 

^ Telescope pointing direction (not necessarily source position) 

^ For these two observations, radio pulsations were detected at a flux density of 20 /xJy (see text). 
Observed 14 times for 1 hour each between this date and 2011-07-11. 



GAMMAY-RAY PULSARS VIA NEW SEARCH METHOD 



TABLE 7 

X-RAY Coverage of the Discovered Gamma-Ray Pulsars 



Pulsar Name 


Instrument 


Exposure Time 


Absorbing column^ 


Flux Upper Limit^ 








(ks) 


(10^1 cm-2) 


(10-13 ergcm-2 s'^) 




J0106+4855 


Suzaku 


23.0 


1.0 


0.843 


>285 


J0622+3749 


Swift XRT 


4.4 


1.0 


2.58 


>65 


J 1620-4927 


XMM-Newton 


6.0 


4.0 


0.674 


>1988 


J 1746-3239 


Swift XRT 


8.7 


1.0 


1.74 


>451 


J 1803-2 149 


Swift XRT 


7.7 


5.0 


3.26 


>34 


J2028+3332 


Swift XRT 


10.3 


1.0 


1.57 


>387 


J2030+4415 


Swift XRT 


10.2 


4.0 


2.53 


>279 


J21 11+4606 


Swift XRT 


10.1 


3.0 


2.25 


>183 


J2139+4716 


Swift XRT 


3.2 


1.0 


3.20 


>78 



Note. — The sky locations of all nine pulsars have been searched for (non-pulsating) X-rays, 
using both archival and new data. No X-ray sources were found at the new pulsar locations, so 
flux upper limits and lower limits on the gamma-ray to X-ray luminosity ratio are reported. 
^ Estimated analogously to Marelli et al.| ( |2011) ). 

^ Upper limit on the unabsorbed flux in the 0.3 - 10 keV energy range, using an absorbed power- 
law model with a photon index of 2 and a signal-to-noise of 3. 
^ Lower limit on the Lj/Lx flux ratio. 
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